Deviations from Full Aging in Numerical Spin Glass Models 
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The deviations from full or pure aging behavior, i.e. perfect t/t w scaling of the correlation and 
response functions of aging glassy systems, are not well understood theoretically. Recent experiments 
of Rodriguez et al. (Phys. Rev. Lett. 91, 037203 (2003)) have shown that full aging applies to the 
thermoremanent magnetization in the limit of infinite cooling rate during the initial quench. In 
numerical models, instantaneous thermal quenches can be — and usually are — applied. This has 
motivated the present numerical investigations of the aging behavior of Edwards- Anderson Ising 
spin glass models with both short- and long-range Gaussian interactions, respectively. We sample 
the distribution of residence time t in suitably defined metastable valleys entered at age t w , finding 
that the deviations from t/t w scaling are small and decrease systematically as the system size grows 
and/or the temperature decreases. Finally, the connection between this behavior and the scaling of 
the correlation function itself is discussed. 

PACS numbers: 02.70.Uu ; 05.40.-a ; 75.50.Lk 
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I. INTRODUCTION 

In a typical aging experiment a spin glass sample is 
quenched below the glass transition temperature in a 
small magnetic field. At time t w after the quench, the 
field is switched off and the subsequent decay of the ther- 
moremanent magnetization is measured. It is well known 
that linear response functions in spin glasses strongly de- 
pend on the age or waiting time t w . In general, the rate 
of change of macroscopic averages in glassy systems sys- 
tematically decreases with t w . 

Considering the broad variety of aging systems, mi- 
croscopic time scales are hardly significant, which high- 
lights t/t w as the obvious choice of scaling variable for 
aging data. The data collapses obtained in this fash- 
ion are fairly good but never perfect: In spin-glass lin- 
ear response experiments P, |3| much better scalings 
are obtained using t/tfo, with the exponent /i slightly 
smaller than one. This situation is referred to as sub- 
aging, since the apparent age tf^ is shorter than t w , as 
opposed to 'pure' or 'full' t/t w aging. In gels || and col- 
loidal glasses jj| , the distribution of the waiting time be- 
tween pairs of consecutive 'anomalous fluctuations' fea- 
tures the opposite deviation, i.e. super- or hyper-aging, 
where \x > 1. Both sub- and super- aging pose a the- 
oretical challenge as they require the introduction and 
interpretation of an additional time scale. This remains 
true even when [i is very close to unity, as in the present 
study. 

Qualitatively, aging reflects the trapping of trajecto- 
ries in metastable 'valleys' of ever increasing thermal 
stability 0, 0, : the pseudo-equilibrium fluctuations 
occur for t < t w within the valley selected at time t w , 
while the exploration of new valleys occurs for t > t w 
and produces the non-equilibrium part of the aging dy- 
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namics. In this picture, the spin rearrangements leading 
from one valley to the next, 'quakes' in our terminology, 
would appear as anomalous events in a fluctuating back- 
ground. The waiting time between quakes corresponds 
to the valley residence time, whose statistical properties 
have repeatedly been considered in models of complex 
relaxation [f| H, ||. According to a recent description 
of non-equilibrium dynamics as a log-Poisson process , 
the distribution R{t | t w ) of residence time t in valleys en- 
tered at time t w should scale as t/t w and possess a finite 
average residence time proportional to t w . The best ob- 
tainable data collapse of the corresponding simulational 
data reveals however a slight super-aging behavior. 

On the experimental side, it has recently been 
shown [T(j that sub-aging of the thermoremanent mag- 
netization can be strongly suppressed by increasing the 
cooling rate of the quench. For the instantaneous quench 
implied in isothermal numerical simulations with a ran- 
dom initial configuration, this strengthens the expecta- 
tion that full aging applies, a conclusion which is not in- 
equivocally supported by the existing numerical studies 
of the correlation function [III Il2| . 

The paper is organized as follows: In Section ^ we 
summarize the exploration method used to identify the 
metastable valleys and describe the procedures followed. 
In Section Hm we present our numerical findings for the 
scaling of the residence time distribution in these valleys 
for spin glass models defined on Euclidean lattices with 
nearest neighbor interactions as well as on random reg- 
ular graphs. In particular, we focus on the system size 
and temperature dependence of [i. Section llVl compares 
the empirical properties of this distribution with the pre- 
diction of the idealized log-Poisson theory of Ref . |7j • We 
show that the theoretical description becomes gradually 
better in the sense that full aging is slowly approached as 
the temperature decreases or the system size increases. 
Finally, the implications for the correlation function are 
discussed, and tentative conclusions are drawn regarding 
the geometrical origin of the deviations from pure aging 
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behavior. 



II. METHOD AND MODELS 

While the fluctuation dynamics of glassy systems 
within metastable valleys is uncontroversially determined 
by the local free energy, a similarly well established theo- 
retical framework for the irreversible quakes interrupting 
the pseudo-equilibrium regime is lacking. Surprisingly, 
there has not been much focus on how these jumps or 
quakes should be identified and characterized. This ques- 
tion is subtle since the first occurrence of a certain rear- 
rangement could well have a non-equilibrium character, 
while subsequent rearrangements of similar nature might 
represent pseudo-equilibrium fluctuations within a larger 
valley. 

The classic analysis of energy landscapes ^3, Ei| sam- 
ples local energy minima using e.g. downhill search al- 
gorithms as thermal quenches, and is not overly con- 
cerned with how attractors are dynamically selected in 
the course of unperturbed isothermal aging. An empir- 
ical procedure to identify the quakes [l5j distinguishes 
'first' and 'subsequent' occurrences of certain events by 
keeping track of the sequence of energies visited while 
aging. Within this sequence, states of energy lower than 
all previously visited states, together with barriers [2(j 
higher than all previously surmounted barriers, define a 
quake as a rearrangement which overcomes a barrier of 
record height and leads to a state of energy lower than 
the lowest energy seen so far. A valley is then simply 
the set of states visited between two subsequent quakes, 
where onset of the latter is marked by the barrier records. 
Very restrictive by construction, this procedure yields no 
output if the dynamics is time translationally invariant, 
and is therefore not available to analyze e.g. the diffu- 
sive regime of supercooled glass formers. Within the ag- 
ing regime, it produces a highly non-trivial yet relati vely 
simple picture of the energy landscape of spin glasses [l6j . 

The Edwards-Anderson model considered below 
is an archetypal glassy system which possesses a com- 
plex dynamics while remaining relatively easy to simu- 
late. The energy of a configuration in this model is given 
by 

E({si,...,s N }) = --^JijSiSj. (1) 

id 

The couplings Jy are symmetric, independent Gaussian 
variables of unit variance. We used standard cubic lat- 
tices with periodic boundary conditions, where Jy is non- 
zero if and only if i and j are nearest neighbors. To in- 
vestigate more general networks as well, we additionally 
explored fc-regular random graphs, i.e. where each spin 
interacts with exactly k other spins chosen at random. 
This is similar to the Viana-Bray model , except that 
the number of links emanating from each spin was fixed 
to k = 6 for better comparison with the Euclidean 3d 
case. 



The systems considered contain up to N — 40 3 spins 
in 3d, and N = 16000 spins in the regular random graph 
case. The dynamics is simulated using the Waiting Time 
Method (WTM) 0. This rejectionless method is equiv- 
alent to the Metropolis algorithm, but much faster at low 
T for simulations of glassy systems. The intrinsic, size in- 
dependent time variable t of the WTM corresponds to the 
number of Monte Carlo (lattice) sweeps in the Metropo- 
lis algorithm as well as to the physical time of a real 
experiment. 



III. NUMERICAL RESULTS 

Using the valley-identification method outlined in Sec- 
tion^ we examine the distribution of residence time t res 
in a valley entered at t w , i.e. 

R(t \t w )=Frob(t res <t\t w ). (2) 

The distribution for a 3c? system with N = 16 3 spins at 
temperature T = 0.5 is shown in Fig^for a wide range of 
t w values. To collapse the data shown we used the scaling 
variable t/tfa with /.i = 1.077(6). The insert confirms that 
full aging, ji = 1, is only a fair approximation. The data 
collapse thus obtained is closely matched by the graph of 
the function f(x) = 1 — (1 + x)~ a , where a ~ 2. The 
latter is not shown in order to avoid masking the data — 
see Fig. 3 of Ref. Q for a direct comparison. 

The present simulations greatly extend our earlier find- 
ings 0, Il6| by probing the size and temperature depen- 
dence of the small deviations from unity of the scaling 
exponent fi. Numerical investigation of such effects re- 
quires the removal of any visible statistical flutter. This 
has been achieved by averaging over at least 4,000 runs 
with different realizations of the couplings Jy for every 
combination of N and T. The calculations were per- 
formed on a large Linux cluster, with a total run-time 
equivalent to four years of computation on a single 2GHz 
processor. The result of these simulations is summed up 
in Fig. 

The high quality of the data collapse in Fig.QJis repre- 
sentative for all sets of runs at fixed size and temperature 
on which the data points in Fig. [21 are based. This sug- 
gests that the 'best' value of n can readily be estimated 
by eye, as small variations of /i immediately produce a 
visible scatter. To remove any personal bias, we did how- 
ever utilize a quantitative error measure e(/i) gauging the 
deviation from perfect collapse of a set of scaled exit time 
distributions: The horizontal standard deviation for a set 
of R(t | t w ) values chosen equidistantly between 0.2 and 
0.8 is calculated for ten values of t/tf^ at fixed fj,. The 
value of e(^) is the sum of these standard deviations, and 
it attains a minimum within a 'reasonable' ^-interval. 
The uncertainty bracket around the [i value producing 
the least scatter is arbitrarily defined by the two values 
of \i for which e(fj.) is no more than 20% larger then its 
minimum value. 
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FIG. 1: The probability distribution of residing at most a 
time t in a valley entered at time t w is calculated for a set of 
10 waiting times equally spaced on a logarithmic scale in the 
interval 10 2 < t w < 10 5 . Each distribution function is plotted 
versus the scaled variable t/t!^, with fi = 1.077, producing 
the perfect data collapse shown in the main panel. All data 
pertain to a 3d spin glass as defined in Section [n] The insert 
shows that choosing /i = 1 gives a visible spread of the curves, 
the latter being shifted to the right as t w grows. The apparent 
age t& inferred from the scaling plots is thus larger than the 
actual age, which is super-aging behavior. 
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FIG. 2: The scaling parameter fj, as a function of temperature 
(left panel) and system size (right panel). For each data point, 
the value of fi is estimated as in Fig. The figure shows that 
super-aging (fi > 1) fades away for low temperatures and 
large system sizes, indicating that deviations from full aging 
stem primarily from finite size effects. 



for large t w values and ruins the scaling. Lacking precise 
estimates of /.i for N < 1000, we present data for larger 
systems only. 



Figure [21 provides information on /i = fj,(N,T) in both 
3d lattices (•) and random graphs (□). The data points 
and error bars are determined as explained above, ensur- 
ing that the quality of the data collapse behind the esti- 
mate of /x is near-perfect, i.e. as good as in Fig.^ In the 
left panel, n(N, T) is shown as a function of T for fixed 
size N = 16 3 . There is a clear decreasing trend, which 
however seems to taper off as the temperature decreases. 
The diminishing deviation from pure scaling concurs with 
the results by Kisker et al. regarding the correlation 
function C (t w ,t w +t) in 3d Gaussian spin glasses. The 
right panel illustrates the behavior of n(N, T) for T = 0.5 
as a function of the logarithm of N . Again, as the number 
of degrees of freedom grows, super-aging is suppressed, 
albeit in a very slow fashion. Finally, we note the striking 
similarity between the Euclidean and random networks 
data in Fig. [21 The scaling of the residence time distribu- 
tion seems to be indifferent to the topology of the graph 
considered. 

When the system size becomes too small, the quality of 
the scaling of R(t \ t w ) begins to deteriorate. This applies 
to systems defined on 3d lattices as well as on random 
graphs. A possible explanation is a finite time (or size) 
effect: by construction, the valley containing the ground 
state will never be exited. By the same token, very low 
lying valleys, which are most likely entered for large t w 
and in small systems, are never exited during the simu- 
lation time. This flattens the s-shape of the distribution 



IV. DISCUSSION 

The residence time analyzed in the previous section 
appears from the outset conceptually similar to the 
trapping time introduced in a heuristic trap model by 
Bouchaud |(| and later considered by Rinn et al. in a 
development specifically dealing with sub- aging behavior. 
However, the statistical properties hypothesized in the 
above models deviate considerably from those presently 
found. Crucially, the fitting form 

R (t i tw ) = i-(i+ t/cr Q (3) 

with a = 2 yields an average residence time 

{tres) = I ll tdt = ^i =t ™> (4) 

which is finite, in contrast with the assumption of the 
trap model [6(- Secondly, the only temperature depen- 
dence is the weak increase of n, while the exponent a 
remains independent of T. 

The numerical results of Fig. [21 points to pure aging 
as the correct limit for small T or large N, although the 
sub-logarithmic size dependence of N prevents us from 
turning this into a firm conclusion. The possible physical 
origin of the deviations from fj, = 1 for finite N was al- 
ready hinted to: Eventually, the trajectory enters the val- 
ley containing the ground state and never leaves. While 
this happens at extreme time scales for realistic system 
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sizes, fewer valleys will in general be available as the tra- 
jectories approach the ground state. The residence time 
for valleys entered at late stages will hence tend to be 
longer than otherwise expected. In accord with this in- 
terpretation, lowering T has an effect on \l similar to the 
effect of enlarging N. In any case, the deviation from 
a pure log-Poisson description of quake dynamics [jj be- 
comes very small for large N, irrespective of the topology 
of the system. This supports the generality of log-Poisson 
statistics, whose prediction for the distribution of the res- 
idence time in a valley entered at age t w coincides with 
Eq. © for n = 1 0. 

While R(t | t w ) might be directly accessible in some ex- 
perimental situations, its relationship to correlation func- 
tions for spin glasses is indirect and requires further as- 
sumptions and/or empirical input. If one assumes with 
Ref. full decorrelation once a trap is left, and that the 
age t w uniquely determines the properties of the 'initial' 
trap, the correlation function basically coincides with the 
probability that no events occur from t w to t w + t. As 
the latter is again given by Eq. ©, this would allow us 
to extend the conclusions already drawn for R to the 
correlation function. 

To assess the validity of such reasoning we will use a 
cartoon picture of the quake dynamics which neglects the 
internal structure of the valleys. This is a questionable 
step, as the time spent searching for the bottom state 
within a valley is of the same order of magnitude as the 
residence time itself [T(| . i.e. quakes are not instantaneous 
events. With this caveat in mind, we write the non- 
equilibrium thermal correlation function as an average 
over the distribution of quakes 0: 

C(t w ,t w +t) = }^P m (t w )P k (t w ,t w +t)c(m,m+k), (5) 

valid in the t > t w regime. In the above formula, P m {t w ) 
is the probability that m quakes are recorded during the 
aging time t w , (t w ,t w + t) is the probability that k 



quakes occurred between t w and t + t w , and, finally, 
c(m, m + k) is the configurational overlap between the 
'bottom states' of the m'th and (m + fc)'th valley. 

The log-Poisson statistics prescribes the form of both 
P m (t w ) and Pk (t w ,t w + t) 7J . While the latter scales 
with t/t w , the former depends on t w alone. This ruins 
pure aging for C(t w ,t w + t) unless the sum of P m {t w ) 
over m happens to factor out. This can only occur if 
the overlap c(m, m + k) does not depend on m. Perfect 
invariance of c(m, m + k) with respect to the index m 
must fail at some stage, since, again, the quake-induced 
rearrangements are likely to change as the ground state 
is approached. Numerical investigations presented else- 
where |7lll6| show that c(m, m + k) has a small but sys- 
tematic m dependence already for small m values. This 
excludes the possibility of 'perfect' pure aging of the cor- 
relation function for the numerical models studied, even 
though pure aging is well fulfilled for the waiting time 
distribution itself. 

As the experimental results of Rodriguez et al. 01 
clearly support pure aging of the correlation function, 
it is possible that a description more refined than what 
Eq. JSJ) provides would modify our conclusion. In view 
of the smallness of the effects involved, it is also possible 
that some subtle dynamical feature of the experimental 
systems is missed by the computer models in the size 
range considered. A final possibility, which we regard as 
the most likely, is that c(m, m + k) is indeed indepen- 
dent of m for systems in the macroscopic limit. This 
would emphasize the extreme non-equilibrium nature of 
the observed low T dynamics in spin glasses. 
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